Mem. S.A.It. Vol. ?, 1 
©SAIt 2008 



Memorie deiia 



1 



2 



Stellar populations 
in tidally stirred dwarf galaxies 

E. L. Lokas', K. Kowalczyk^ and S. Kazantzidis^ 

Nicolaus Copernicus Astronomical Center, Bartycka 18, 00-716 Warsaw, Poland, e-mail: 
lokas@caink. edu.pl 

Warsaw University Observatory, Al. Ujazdowskie 4, 00-478 Warsaw, Poland 

Center for Cosmology and Astro-Particle Physics and Department of Physics and 

Department of Astronomy, The Ohio State University, Columbus, OH 43210, USA 



Abstract. Using A'-body simulations we study the evolution of separate stellar populations 
in dwarf galaxies in the context of the tidal stirring scenario for the formation of dwarf 
spheroidal (dSph) galaxies in the Local Group. The dwarf galaxies, initially composed of 
a stellar disk and a dark matter halo, are placed on seven different orbits around the Milky 
(~ i . Way. The stars are divided into two populations, within and outside the half-light radius, and 

^ \ their positions are followed for 10 Gyr. We find that the populations retain different density 

distributions even over such long timescales. Some of the stars of the outer population 
migrate to the central part of the dwarf forming an extended core while the stars of the inner 
population develop a tail in the outer parts. In addition, the outer population is more heavily 
stripped by tidal forces from the Milky Way and may become subdominant at all radii on 
tight enough orbits. We conclude that the tidal stirring model is fully compatible with the 
presence of multiple stellar populations in dSph galaxies. 
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1. Introduction tified were shown to be distributed differently 

in space and possess various velocity distribu- 
Dwarf spheroidal (dSph) galaxies of the Local ^^^^^ ^s quantified e.g. by their velocity disper- 
Group provide unique laboratory for studying ^^^^ profiles (Tolstoy et al. 2004; Koch et al. 
the processes of galaxy formation and evo- 2OO6; Battaglia et al. 2008; Kirby et al. 2011; 
lution. Their relative proximity has recently je Boer et al. 2012). The presence of multiple 
enabled detailed probing of density distribu- stellar populations has been invoked as a possi- 
tions and kinematics, as well as metallicities of ble explanation for anomalous velocity disper- 
their stellar components. These observational gion profiles observed in some dSph galaxies 
studies have revealed the presence of multiple (j^ata et al. 2006; McConnachie et al. 2007). It 
stellar populations in the dwarfs, differing by has also recently proved useful in constraining 
metallicity and age. The populations thus iden- the density slope of dSph galaxies (Walker & 
Penarrubia 20 11 ) . 
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Table 1. Parameters of the simulated dwarfs 
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Substantial progress has also been achieved 
in theoretical modelling of the formation of 
dSph galaxies. One promising scenario for 
their origin invokes tidal evolution in the vicin- 
ity of more massive galaxies like the Milky 
Way or Andromeda (Mayer et al. 2001). The 
tidal stirring model envisions the present- 
day dSph galaxies as products of the re- 
peated action of the tidal force from the host 
galaxy on a late-type progenitor accreted a 
few Gyr ago, thus explaining naturally the 
observed morphology-density relation in the 
Local Group. The process can be rehably mod- 
elled using A^-body simulations which demon- 
strate that such evolution is indeed feasible 
over reasonable timescales (Klimentowski et 
al. 2009; Kazantzidis et al. 2011). The transfor- 
mation from disks to spheroids is accompanied 
by randomization of stellar orbits and strong 
mass loss due to tidal stripping. 

The morphological evolution typically in- 
volves a bar-hke stage (Lokas et al. 2012) and 
therefore a question arises what effect the in- 
stability associated with bar formation may 
have on the properties of stellar populations 
present in the dwarf galaxy. Since most of the 
stars in dSph galaxies are old, they have prob- 
ably formed before the dwarf's progenitor was 
accreted by the host. Although some star for- 
mation may occur as a result of tidal stirring, 
the amount of new stars is probably not large 
as the gas reservoir is likely to be depleted 
quickly by ram pressure stripping (Mayer et 
al. 2007). In this work we address the ques- 



tion whether the separate stellar populations 
present in the progenitor prior to accretion are 
likely to survive until the present time. Such 
considerations may put some constraints on the 
amount of tidal stirring experienced by dSph 
galaxies in our vicinity. 



2. The simulations 

For the present study we used a subset of 
tidal stirring simulations described in detail 
in Kazantzidis et al. (2011) and Lokas et al. 
(2011). Our dwarf galaxy was initially com- 
posed of an exponential stellar disk of mass 
2 X IO^Mq embedded in a dark matter halo 
of virial mass 10^ M0 and an NFW (Navarro 
et al. 1996) density profile with concentration 
c = 20. The A'^-body realization of the dwarf 
contained 1.2 x 10'' stars and 10^ dark mat- 
ter particles. The dwarf galaxy was placed on 
seven different orbits (initially at the apocen- 
ter) around a host galaxy with the present- 
day properties of the Milky Way (Widrow & 
Dubinski 2005) and its evolution was followed 
for 10 Gyr with the PKDGRAV code (Stadel 
2001). 

The orbital parameters of the simulations 
are listed in Table 1 . Columns 2-6 of the Table 
give respectively the apocenter r-^^o, the peri- 
center Tperi, the ratio of the two distances, the 
orbital time and the number of pericenters 
passed during the 10 Gyr of evolution. The or- 
bits cover a large range of sizes and eccen- 
tricities resulting in different amounts of tidal 
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Fig. 1. Evolution of the density profiles of stars in the inner Population 1 (dashed lines) and the outer 
Population 2 (solid lines) for the dwarf galaxy on orbit Ol. The populations arc initially divided by the half- 
light radius (upper left panel). The panels (from the upper left to the lower right) show the density profiles 
at subsequent apocenters of the orbit and at the final stage. The times of the measurements are given in each 
panel. 



force experienced by the dwarf galaxy. The last 
three columns of Table 1 list the observational 
parameters of the dwarfs at the final stage as 
they would be measured by an imaginary ob- 
server at the center of the Milky Way: the abso- 
lute magnitude, the projected half-light radius 
and the central surface brightness (Lokas et al. 
201 1). The parameters fall in the ranges occu- 
pied by the classical dSph galaxies in the vicin- 
ity of the Milky Way. 

3. Evolution of stellar populations 

In order to see how stellar populations may be 
affected by tidal evolution we divided the stel- 
lar component of the dwarf galaxy into two 
populations, each containing initially half of 
the total number of stars. The inner population 
(Population 1) was selected to contain all stars 
with radii r < r\j2 and the outer population 



(Population 2) all stars with r > ri/2 where 
ri/2 = 0.7 kpc is the three-dimensional half- 
light radius of the initial stellar disk. Although 
quite simpUstic, the division probably reflects 
to some extent the actual distribution of stars 
of different metallicity following from multiple 
star formation episodes. The stars formed later 
are likely to be produced from the gas expelled 
at earlier stages which was later accreted and 
sunk to the center of the dwarf galaxy resulting 
in the segregation of populations with different 
stellar properties. 

The positions of the stars in each popula- 
tion were traced for the whole evolution time. 
Figure 1 shows the density profiles of the stars 
in the two populations at the initial time, the 
subsequent apocenters and at the final out- 
put for orbit 01, which may be considered a 
typical orbit of dwarf galaxies in the vicin- 
ity of the Milky Way (Diemand et al. 2007; 
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Fig. 2. The density profiles of stars in the inner Population 1 (dashed lines) and the outer Population 2 
(solid lines) at the end of evolution (after 10 Gyr) for orbits 02-07 (from the upper left to the lower right 
panel). The populations were initially divided by the half-light radius (as in the upper left panel of Figure 1). 



Klimentowski et al. 2010). Although the stars 
of the two populations were initially spatially 
separated, as the evolution proceeds the outer 
population tends to populate the inner region 
forming a core-Uke distribution, while the in- 
ner population extends to larger radii forming 
a tail. The stars are lost from both popula- 
tions due to tidal forces, as confirmed by the 
decreasing normalization of the density pro- 
files. For most of the time. Population 1 dom- 
inates in the center, and Population 2 in the 
outer parts. However, with time the outer pop- 
ulation is more heavily stripped and gradually 
becomes subdominant at all radii. In addition, 
both components display a flattening of their 
respective density profiles at radii r > 1 kpc 
signifying the transition to tidal tails. 

As expected, the final outcome of the evo- 
lution depends sensitively on the orbit on 
which the dwarf is accreted. We illustrate this 
in Figure 2 which shows the final density pro- 



files of the two populations on orbits 02-07. 
The populations were defined initially in ex- 
actly the same way as for orbit 01. We see 
that on the very tight orbit 02 and the one with 
the smallest pericenter 04 the outer population 
is very strongly stripped and very few stars of 
this population are still bound to the dwarf. On 
the other hand, the less tight orbits 03, 05 and 
06 allow both populations to survive and retain 
their different distributions. The most eccentric 
orbit 07 is an interesting intermediate case. 

4. Conclusions 

Using A^-body simulations of dwarf galax- 
ies orbiting the Milky Way we have studied 
the evolution of separate stellar populations in 
dwarfs. We have shown that in spite of re- 
peated action of tidal forces at pericenters, the 
differences in the distribution of stellar pop- 
ulations survive over timescales comparable 
to the present age of the Universe. The stars 
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of the inner stellar component migrate out- 
wards and the ones of the outer one inwards 
populating the regions initially devoid of the 
stars of each type. Still, the populations do not 
mix completely and the density profiles of the 
two components remain distinguishable. The 
inner component remains more centrally con- 
centrated and the outer one more extended so 
that their respective characteristic radii (e.g. 
the half-Ught radii) at present may be very dif- 
ferent. In general, the more extended popula- 
tion is stripped more effectively and thus may 
become subdominant at all radii for sufficiently 
tight orbits. 

We have demonstrated that the presence 
of multiple stellar populations in dSph galax- 
ies of the Local Group is entirely consistent 
with the tidal stirring scenario which envisions 
their formation as following from interactions 
of late-type progenitors with a Milky Way-like 
host galaxies. Even with the instability inher- 
ent to the bar formation usually taking place at 
the first pericenter passage, the separate pop- 
ulations do survive. This should be the case 
also for the generic example of such inter- 
actions, namely the Sagittarius dwarf, which 
moves around the Milky Way on a rather 
tight orbit and is thus subject to particularly 
strong tidal effects resulting in the formation of 
pronounced tidal tails (Majewski et al. 2003; 
Law et al. 2005; Lokas et al. 2010). Indeed, 
the presence of multiple stellar populations in 
Sagittarius and its tidal tails is well established 
observationally (Siegel et al. 2007; Chou et al. 
2007). 

Acknowledgements. This research and ELL's par- 
ticipation in EWASS 2012 was partially supported 
by the Polish National Science Centre under grant 
NN203580940. KK acknowledges the summer stu- 
dent program at the Copernicus Center in Warsaw. 
We are grateful to A. Helmi for the suggestion to 
investigate the evolution of separate stellar popula- 
tions in tidally stirred dwarfs. 

References 

BattagUa, G., Hehni, A., Tolstoy, E., et al. 

2008, ApJ, 681, LI 3 
Chou, M.-Y., Majewski, S. R., Cunha, K., et 

al., 2007, ApJ, 670, 346 



de Boer, T. J. L., Tolstoy, E., Hill, V., et al. 

2012, A&A, 544, A73 
Diemand, J., Kuhlen, M., & Madau, P. 2007, 

ApJ, 667, 859 
Ibata, R., Chapman, S., Irwin, M., Lewis, G., 

& Martin, N. 2006, MNRAS, 373, L70 
Kazantzidis, S., Lokas, E. L., Callegari, S., 

Mayer, L., & Moustakas, L. A. 2011, ApJ, 

726, 98 

Kirby, E. N., Lanfranchi, G. A., Simon, J. D., 
Cohen, J. G., & Guhathakurta, P. 2011, ApJ, 

727, 78 

Klimentowski, J., Lokas, E. L., Kazantzidis, S., 
Mayer, L., & Mamon, G. A. 2009, MNRAS, 
397, 2015 

Klimentowski, J., Lokas, E. L., Knebe, A., et 

al. 2010, MNRAS, 402, 1899 
Koch, A., Grebel, E. K., Wyse, R. R G., et al. 

2006, AJ, 131,895 
Law, D. R., Johnston, K. V., & Majewski, S. R. 

2005, ApJ, 619, 807 
Lokas, E. L., Kazantzidis, S., Majewski, S. R., 

etal. 2010, ApJ, 725, 1516 
Lokas, E. L., Kazantzidis, S., & Mayer, L. 

2011, ApJ, 739, 46 
Lokas, E. L., Majewski, S. R., Kazantzidis, S., 

etal. 2012, ApJ, 751,61 
Majewski, S. R., Skrutskie, M. P., Weinberg, 

M. D., & Ostheimer, J. C. 2003, ApJ, 599, 

1082 

Mayer, L., Govemato P., Colpi, M., et al. 2001, 

ApJ, 559, 754 
Mayer, L., Kazantzidis, S., Mastropietro, C, & 

Wadsley, J. 2007, Nature, 445, 738 
McConnachie, A. W, Penarrubia, J., & 

Navarro, J. R 2007, MNRAS, 380, L75 
Navarro, J. R, Frenk, C. S., & White, S. D. M. 

1996, ApJ, 462, 563 
Siegel, M. H., Dotter, A., Majewski, S. R., et 

al. 2007, ApJ, 667, L57 
Stadel, J. G. 2001, PhD thesis, Univ. of 

Washington 
Tolstoy, E., Irwin, M. J., Hehni, A., et al. 2004, 

ApJ, 617, LI 19 
Walker, M. G., & Peiiarrubia, J. 2011, ApJ, 

742, 20 

Widrow, L. M., & Dubinski, J. 2005, ApJ, 631, 
838 



